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Abstract 

The tensor polarization of particles and nuclei becomes constant in the coordinate system rotat¬ 
ing with the same angular velocity as the spin and rotates in the lab frame with the above angular 
velocity. The general equation defining the time dependence of the tensor polarization is derived. 
An explicit form of dynamics of this polarization is found in the case when the initial polarization 
is axially symmetric. 
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I. INTRODUCTION 


Polarized beams are often considered as a research tool to study the fundamental interac¬ 
tions and to search for new physics. Unlike spin-1/2 fermions, particles and nuclei with spin 
s > 1 possess tensor polarization. This polarization is an important property of particles 
and nuclei and can be measured with a good accuracy. Investigations of dynamics of ten¬ 


sor polarization o 
[l, y] , to predict 


it nuclei (e.g., the deuteron with spin 1) made it possible to discover 


4j, and to describe [5H7J some new effects. However, these effects are 


conditioned either by quadratic in spin interactions of deuteron 


-6] or by quantum beats 


in positronium [7]. To get a basic understanding of polarization effects, it is necessary to 
consider linear in spin interactions of tensor-polarized, particles and nuclei. While this prob¬ 
lem is rather important, it is not much discussed. For all we know, the dynamics of tensor 
polarization in external fields has been considered in Ref. jg]. In this work, some examples 
of evolution of the tensor polarization have been investigated and some properties of the 
5x5 spin transfer matrix describing this polarization have been considered. Evidently, the 
use of the 5x5 matrix is much more difficult than that of the corresponding 2x2 matrix 
for the vector polarization. 

In the present work, we develop a simple general approach to a description of dynamics 
of the tensor polarization of particles and nuclei when this dynamics is caused by linear in 
spin interactions with external fields. The approach proposed allows us to couple dynamics 
of the tensor and vector polarizations and to derive general formulas defining an evolution 
of particles/nuclei polarization in external fields. 


II. GENERAL PROPERTIES OF SPIN DYNAMICS 


As is known, the spin rotation is exhaustively described with the polarization vector P 
defined by 


p < ^ > . 

Pi = -, i = x,y,z. 


( 1 ) 


Here Si are corresponding spin matrices and s is the spin quantum number. Averages of spin 
operators are expressed by their convolutions with the wave function T(f). Particles with 
spin s > 1 also possess a tensor polarization. Main characteristics of such a polarization are 
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specified by the polarization tensor P l3 , which is given by 


3 < SiSj + SjSi > —2 s(s + 1 )Sij 
Pn = --——T\--> hJ=x,y,z. 


( 2 ) 


2s(2s - 1) 

The polarization tensor satisfies the conditions Pij = Pji and P xx + P yy + P zz = 0 and 
therefore has five independent components. In the general case, the polarization vector and 
the polarization tensor are time-dependent. Additional tensors composed of products of 
three or more spin matrices are needed only for the exhaustive description of polarization 
of particles/nuclei with spin s > 3/2. 

The eight parameters defined by the polarization vector and the polarization tensor are 
independent. In particular, < SiSj >^< Si >< S 3 >. 

The quantum-mechanical Hamiltonian describing the interaction of the spin with external 
fields contains terms linear and bilinear on the spin: 


n = n-s + Q jk s j s k . ( 3 ) 

Here S is the matrix operator describing the rest-frame spin. Since the spin matrices satisfy 
the commutation relation 


[Si, S/] = ietjkSk ( i,j , k = x,y, z ), (4) 

the first term in Eq. (J3J) defines the spin rotation with the angular velocity f2. The absolute 
value and the direction of may arbitrarily depend on time. 

Even when the second term in Eq. (J3j) is omitted, dynamics of the polarization tensor is 
not trivial. We cannot characterize the particle spin wave function only by a single frequency. 
For a spin-1 particle, there are two frequencies, and —H, defining three equidistant levels. 
For a particle with the integer spin s = N, there are 2N frequencies, ±f2, ±2f2, ..., ±NQ, 
defining 2N + 1 equidistant levels. The single frequency of the spin rotation, 0, originates 
from the specific commutation relations (J4]) for the spin operators S For a spin-1 particle, 
this property of the spin operators manifests in the fact that (Si)i 3 = (<S'j) 3 i = 0 (i = x, y, z). 
Therefore, the spin operators mix only two neighboring levels. But it is not the case for the 
tensor polarization operators. Three of them mix the levels 1 and 3 because the components 
(^)i 3 , ( 5 ^) 31 , ( 5 ^) 13 , ( 5 ^) 31 , (. S x Sy + SyS x ) 13 , and (. S x S y + S y S x ) 31 are nonzero. As a 
result, the frequency 2also appears. For particles with greater spins (s > 1), the situation 
is still more complicated. 
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We consider electromagnetic interactions, namely, coupling of spins to electric and mag¬ 
netic fields. The second term in Eq. (J3J) characterizes the tensor electric and magnetic 
polarizabilities and the quadrupole interaction. It influences dynamics of both the polar¬ 
ization vector and the polarization tensor. As a rule, this term is much less than the first 
one. Its order of magnitude can be easily determined. We can calculate the commutator 
[H, S'*] taking into account that s ~ 1. The term under consideration brings a correction 
of the order of \QjkSjSk\ to the angular frequency of spin precession hi. Let us evaluate 

on spin dynamics is rather complicate and does not reduce to a 
). The angular frequency of spin precession of the deuteron in 


this correction whose effec 
change of f2 (see Refs. 


the magnetic field B — 3 T is equal to 1.2 x Hr s . If the tensor magnetic polarizability of 
the deuteron is close to the value (3t = 1.95 x 1CT 40 cm 3 predicted by the theory [r], Q, the 
related effective change of the angular frequency of spin precession is about 2 x 1CT 4 s -1 . 
The predicted values of the tensor electric polarizability of the deuteron 15H171] are slightly 


less and cannot lead to a greater effect. On the other hand, spin-tensor effects can be caused 
by intrinsic electric quadrupole moments of nuclei. For example, a spatial derivative of an 
electric field strength in Penning traps is usually of the order of 10 6 V/m 2 . The quadrupole 
moment of the deuteron is equal to Qd = 0.286 fm 2 . A simple estimate shows that the 
related effective change of the angular frequency of spin precession of the deuteron is of 
the order of 1CT 8 s -1 . For polarized deuteron beams in accelerators and storage rings, the 
intrinsic electric quadrupole moment also causes spin-tensor interactions. However, its effect 
on the spin evolution is a few orders of magnitude lower than that of the tensor magnetic 
and electric polarizabilities. When B = 3 T, the ring radius R ~ 10 m, and the held index 
n = —(R/B)(dB/dr) ~ 0.1, the effective change of the angular frequency of spin precession 
caused by the electric quadrupole moment is of the order of 1CT 9 s _1 even for relativistic 
deuterons. 

Thus, one can usually neglect an influence of the spin-tensor interactions on spin dy¬ 
namics. Such an influence can be significant only in special cases like hyperfine interactions 
in atoms or in experiments specially designed to measure the tensor magnetic and electric 
polarizabilities. Hereafter, we will disregard the spin-tensor interactions and will take into 
account only the first term in the interaction Hamiltonian (J3j). 
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The equation of spin motion is given by 


g = i[«,3 = nx5. 


(6) 


When the particle/nucleus moves in electromagnetic fields, f2 is defined by the Tliomas- 
Bargmann-Michel-Telegdi (T-BMT) equation (9j extended on the electric dipole moment 
(EDM): 


n = 


e 

me 


G + 
*1 


7 


B 


7 G 


7 + 1 


(G + 


E- 


7 


7 + 1 


({3-E)f3 + /3xB 


7 + 1 


(3xE 


( 6 ) 


Here G = (g— 2)/2, g = 2 men/ (efts), g = 2 med/ ( ehs ), f3 = v/c, 7 is the Lorentz factor, and 
g and d are the magnetic and electric dipole moments. An extension of the T-BMT equation 
due to the EDM has already been discussed in the original paper of Bargmann, Michel and 


Telegdi 


9]. Then, the equation of spin motion of the partic 


e with the anomalous magnetic 
El| by the dual transformation 


moment // and the EDM has been obtained in Refs. jlQ . 
g! —> d, B —>■ E, E —>■ —B. The rigorous derivation of this equation has been presented 
in Ref. jl2]. The resulting equation of spin motion coincides with that presented in Refs. 
jlQ . 11 ]. However, the derivation fulfilled in Ref. 12] has not used the supplementary 
assumption of dual symmetry. 

The Pauli matrices together with the unit matrix describe the spin of a Dirac particle 
and generate an irreducible representation of the 517(2) group. Algebraically, the 517(2) 
group is a double cover of the three-dimensional rotation group 50(3). As a result, the 
spin dynamics defined by the Dirac equation fully corresponds to the classical picture of 
rotation of an intrinsic angular momentum (spin) in external fields. The angular velocity of 
spin rotation does not depend on the spin quantum number. Therefore, both the classical 
description and the formalism based on the Pauli matrices are applicable to particles/nuclei 
with an arbitrary spin if the effect of spin rotation is analyzed. Certainly, this formalism 
becomes insufficient, if one considers spin-tensor effects mentioned in Sec. U 

In Eq. (| 6 ]), we can estimate terms dependent on the EDM. Any experimental constraints 
on the deuteron EDM are not established yet. For the proton, the current indirect bound 
of \d p \ < 7.9 x 10 ~ 25 e-cm has been obtained using Hg atoms 18]. If the deuteron EDM was 
equal to this value, its contribution to the angular frequency of spin precession of 1 GeV/c 
deuterons in the magnetic field B = 3 T were equal to 5 x 10 -3 s -1 . This means that one 
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can often neglect an influence of the spin-tensor terms in the Hamiltonian (j3J) on the spin 
dynamics but can keep the EDM-dependent terms in Eq. (ED- Certainly, a possible future 
ascertainment of very strong restrictions on the deuteron EDM can change this situation. 


III. HUANG-LEE-RATNER APPROACH 


Let us consider the approach used by Huang, Lee, and Ratner p] to describe an evolution 
of the tensor polarization. Instead of the polarization tensor (J2]), the following spin-tensor 
operators can be used 


To = -L (3SJ - 2). T ±1 =±Y(S±S 3 + S 3 S ± ), T ±2 =LL|. (7) 

To obtain observable quantities, these operators should be averaged. 

The approach proposed in Ref. jsj was based on the T-BMT equation and disregarded the 
spin-tensor interactions in the Hamiltonian (J3J). This approach used the the Frenet-Serret 
curvilinear coordinates. Expressing the spin vector in terms of its components, 


S — Siei + S2C-2 + *S'3 e 3, 


( 8 ) 


and defining S± — S\ ± 1 S 2 , one obtains the following equation of spin motion: 
dS± 


dqb 


= FiG^S± ± iF±Ss, 


f = i< F - s +- F + s ->- 


( 9 ) 


where F± characterize the spin depolarization kick and (f) is the azimuthal orbit rotation 
angle. The spin tune, Gy, is the number of spin revolutions per orbit turn. 

Equations (ElR© allow one to obtain the set of equations defining the evolution of the 
tensor polarization 8]: 


d7o_J_/R?3 0 , 0 dS 3 \ .y/E 


d(f) \/2 \ d(f> 3 F 3 d(j> 


= (F_T + i + F + T U), 


dT- 


±1 


\/3 (dS± 0 , c dS± | dS 3 c c , dSz\ , ^ ^ 

d<j> 2 V d(j> S * + dcj) + d(f) S± + S± d(j) ) ±lGlT±l 1 2 F±T ° ^ 


dT. 


±2 


d(f) 


— f ^ S± + ) = ±2*G 7 T ±2 - iF ± T ±] . 


2 \ dcp 


dcj) 


( 10 ) 
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These equations can be presented in the matrix form: 


dT 

# “ ’ 


T = 


( T +2 ^ 


l 2 iGy 

-iF + 

0 

Vq 

0 

0 ' 

T+i 


-iF _ 

iGy 


0 

0 

T 0 

, 71 = 

0 

~ % ~2 F ~ 

0 

~ l ^2 F+ 

0 

T i 







— ± 


0 

0 


—iGy 

~iF + 

T- 2 j 


l 0 

0 

z 

0 

- iF _ - 

-2iGy j 

some applications 

of this approach have been considered. 



( 11 ) 


In Ref. |s |, 

While the approach of Huang, Lee, and Ratner is absolutely correct, we can make some 
evident remarks. First of all, the use of 5 x 5 matrices instead of standard 2x2 acceler¬ 
ator matrices needs much more cumbersome derivations. Moreover, even investigations of 


spin-tensor effects performed in Ref. Q] used three-component wave functions and 3x3 
matrices. Another remark is an absence of a direct connection between evolutions of the 
polarization vector and the polarization tensor, while spin-tensor interactions are not taken 
into consideration in Ref. js]. 

In the present work, we propose a different approach which couples dynamics of the 
tensor and vector polarizations and simplifies a description of the evolution of the tensor 
polarization of partcles/nuclci in external fields. 


IV. EVOLUTION OF THE POLARIZATION TENSOR 

It is not evident whether the evolution of the polarization tensor also reduces to the 
rotation. To clear the problem, we use the following approach. Let us consider the Cartesian 
coordinate system rotating about the z axis with the same angular velocity Q(t) as the spin. 
In the general case, this angular velocity depends on time. We may superpose the rotating 
and nonrotating coordinate systems at the initial moment of time t — 0 (e)(0) = e*). The 
rotating coordinate system is denoted by primes. 

The spin components in the rotating coordinate system remain unchanged: 

-± = 0 (i = x ,y,z). (12) 
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As a result, all tensor polarization operators and all components of the polarization 
tensor are also unchanged in this coordinate system (S^Sj + SjS^ = const, P'j = const). 
This important property shows that the tensor polarization of particles/nuclei with spin 
s > 1 rotates in external fields similarly to the vector polarization. This is valid not only 
for electromagnetic interaction but also for other (weak, gravitational) interactions. Other 
products of the spin operators S^S) .. . S' k are also conserved in the rotating coordinate 
system. 

To determine the particle polarization in the nonrotating coordinate system (lab frame), 
we need to connect the directions of basic vectors of the primed and unprimed coordinate 
systems. This problem can be easily solved. We choose the initial moment of time t — 0. 
Evolution of any basic vector of the primed coordinate system coincides with that of the spin 
when it is initially directed along the considered basic vector. Therefore, the final directions 
of the three primed basic vectors are defined by the three final directions of the spin when its 
corresponding initial ones are parallel to the three Cartesian axes of the unprimed coordinate 
system. As a result, dynamics of the polarization tensor reduces to that of the polarization 
vector. In particular, the evolution of the polarization tensor in accelerators and storage 
rings can be unambiguously defined with usual 2x2 spin transfer matrices. Thus, one need 
not apply more cumbersome 5x5 matrices. 

Time dependence of the polarization tensor defined in the lab frame can be easily ex¬ 
pressed in terms of the basic vectors e\. Since e'(0) = e, and 

s'(t) = E (<=:«' ( 13 > 

k 

the polarization tensor is given by 

W) = 2.,( 2 L 1) E [(«{(*)' «»)(« 5 (*)' «) < StS, + S,s k >] - ±±iv ( 14 ) 

This simple equation defines general dynamics of the tensor polarization of particles and 
nuclei in external fields. 

V. EVOLUTION OF AXIALLY SYMMETRIC POLARIZATION 

A calculation of the polarization tensor can usually be still more simplified. In particular, 
this is possible when the initial polarization is axially symmetric. It is convenient to direct 



the z' axis along the symmetry axis. For the considered initial polarization, P' xx = P' yy 
and P'j = 0 when i ^ j. It has been ascertained in the precedent section that the primed 
polarization tensor remains unchanged. 

Let us suppose that the direction of e' z (t ) in the lab frame is defined by the time-dependent 
spherical angles 9 and 0. In the considered case, the lab frame polarization tensor does not 
depend on directions of the x' and y' axes in the plane orthogonal to z' and we can choose 
the primed basic vectors as follows: 


e' x = e x cos 9 cos 4> + e v cos 6 sin (f> — e z sin 6 , e' = — e x sin 0 + e v cos 0, 


(15) 


e' z = e x sin 6 cos (j) + e y sin 6 sin 0 + e 2 cos 6. 

The connection between the spin operators in the unprimed and primed coordinate sys¬ 
tems is defined by Eqs. (1T3|) and (TT5|) and is given by 

S x = S' x cos 9 cos 0 — S' y sin 0 + S' z sin 9 cos 0, S y = S' x cos 9 sin 0 + S' y cos 0 + S' z sin 9 sin 0, 

S z = — S' x sin 9 + S' z cos 9. 

(16) 

Due to the axial symmetry, averaging vanishes < S' x > and < S' y >. The polarization vector 
P =< S > /s is given by 

P x (t ) = P 2 (0) sin 9 cos 0, P y (t ) = P 2 (0) sin 0 sin 0, P z (t ) = P 2 (O)cos0. (17) 

The use of Eqs. (111)) and dill) leads to the following equation defining the polarization 
tensor: 


P xx (t) = A ( sin 2 9 cos 2 0 — 


P yy (t ) = A I sin 2 9 sin 2 0 — 


P zz {t ) = A l^cos 2 9 — - J , P xy (t ) = A sin" 9 sin 0 cos 0, 

Pxzit) = A sin 9 cos 9 cos 0, P yz (t) = A sin 9 cos 9 sin 0, A = -P zz (0). 


(18) 


This general equation exhaustively describes dynamics of the tensor polarization of par¬ 
ticles and nuclei in external fields. The time dependence of the angles 9 and 0 can be 
determined with the use of the spin motion equations (E]) and ([6]) . For particles/nuclei in 
accelerators and storage rings, 2x2 spin transfer matrices are commonly applied 0 Q- 
To take into account spin-tensor interactions transforming the tensor polarization into the 
vector one and the other way round, 3x3 spin matrices should be used (see Ref. js|). 
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VI. SUMMARY 


In the present work, we have considered the evolution of the tensor polarization of par¬ 
ticles and nuclei when this evolution is caused by linear in spin interactions with external 
fields. The tensor polarization becomes constant in the coordinate system rotating with the 
same angular velocity as the spin. This means that the tensor polarization rotates in the lab 
frame with the above angular velocity. We have derived the general equation (fl4l) defining 
the time dependence of this polarization. Calculations becomes appreciably simpler when 
the initial polarization is axially symmetric. In this case, the explicit form of dynamics of 
the tensor polarization is given by Eq. (fTHjh 

The approach used has allowed us to couple dynamics of the tensor and vector polariza¬ 
tions. This result considerably simplifies a determination of the tensor polarization because 
the vector one is exhaustively described by 2 x 2 spin transfer matrices. As a result, one need 
not apply 5x5 matrices to investigate the evolution of five independent components of the 
polarization tensor. The obtained general equations can be widely used (e.g., for polarized 
beams in accelerators and storage rings). 
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